INVARIANCE AND INNER FRACTALS IN POLYNOMIAL AND 
TRANSCENDENTAL FRACTALS 
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Abstract. A lot of formal and informal recreational study took place in the 
fields of Meromorphic Maps, since Mandelbrot popularized the map z 4— z 2 + c. 
An immediate generalization of the Mandelbrot z <— z n + c also known as the 
Multibrot family were also studied. In the current paper, general truncated 
polynomial maps of the form z <— ^2 p>2 a p xP + c are studied. Two fundamental 
properties of these polynomial maps are hereby presented. One of them is the 
existence of shape preserving transformations on fractal images, and another one 
is the existence of embedded Multibrot fractals inside a polynomial fractal. Any 
transform expression with transcendental terms also shows embedded Multibrot 
fractals, due to Taylor series expansion possible on the transcendental functions. 
We present a method by which existence of embedded fractals can be predicted. 
A gallery of images is presented alongside to showcase the findings. 



1. Introduction 

Mandelbrot popularized the map T(z) : z <— z 2 +c which bears his name (p],[5]). 
The Mandelbrot set M, comprise of points L z\ such that T(z), repeated infinite 
time, upon itself, remains bounded. 

The iteration is defined as in pQ :- 

(1) Initialize current point c = P(x,y) and z = c. 

(2) Use the Transform map T : z <— z 2 + c to generate z\. 

(3) Repeat the procedure (2), to get z 2 ,z 3 , ...,zi until \z\\ > k or I = N is the 
maximum number of iteration. 

(4) If iterations expired, the original point c G M. 

(5) Else c M. 

(6) Take another point in (x,y) plane, and start again from step (1), till all 
the points in the plane are exhausted. 

When T(z) is applied on itself once it can be said T 2 (z). Making things general, 
an 'x' times repetition of T{z) upon itself would be termed as T x (z). 
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Then, Formally:- 

M(k) = {z : lim \T x (z)\ < k} (1.1) 

x— >oo 

A plot of the set yields outstanding visual imagery, and patterns. 

1.1. Escape Time Fractals. A generic way to produce any set $ of this kind 
would be then to have a transform function F(z) : z F(z) such as to:- 

$(fc) = {z : lim \F x {z)\ < k} (1.2) 

x— >oo 

Clearly then, k is a very important parameter. However the removal of k can 
be easily done with:- 

$ = {z : lim \F x (z)\ < oo} (1.3) 

For in case of Mandelbrot set M it can be easily shown that any point \z\ > 2 
would escape to infinity, that means:- 

M(z 2 + c, 2) = M(z 2 + c,k); 2 < k < oo 

so for M , k is actually omitted. 



The images generated by equation (1.3) are called "Escape Time Fractals", also 
known as "orbits" fractals. 

1.2. Our Methodology. We have modified an open source 'Mandelbrot Set Viewer' 
computer program, and incorporated the ability to type in mathematical expres- 



sions for the function F(z) in equation (1.3). The source code is available here 



There is an expression handler in the program that can handle many standard 
mathematical functions. Given the expression, the program would evaluate the 
expression dynamically, to generate the 'escape time fractal' for the function F(z). 
The program is capable of:- 

(1) zoom-in and zoom-out. 

(2) modifying the radii of convergence l k\ to be expressed as a power of 'e'. 

(3) modifying number of iterations to perform as 'infinity'. 

(4) saving the image as 'png' format image file, that is how illustrations of this 
paper were generated. 

(5) plotting the set in Colors, to be discussed more in next. 

1.2.1. Plotting the Set using Software. From a computational point of view, there 



is no way to perform the iteration of equation (1.3) up-to infinity. Hence, a max- 
imum W number of iteration is performed. Using the custom software we can 
manipulate i N\ N G [0, 99999]. If a point z exceeds the limit k at iteration number 
m, then a ratio x = m/N is defined. Clearly x — 1 means, the function remained 
bounded till the 'iV'th iteration. Based upon the value of x, the color of the point 
is chosen, among a fixed set of colors. Though a point might show x = 1.0, that 
does not imply that, it would remain so for an increased N. Also, another point 
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having current x < might get into the set, for an increased k. The coloring can 
be formally defined as a mapping:- 

C(x) : R ->■ N (1.4) 

The custom application which produced the illustrations of present paper is ca- 
pable of plotting in any integer number of colors, and can handle any mathematical 
expression and function. The color codes are dependent on the color map used. 



The program uses equation (1.5) for coloring 



C{x) : \p x] ; Vx G (0, 1) (1.5) 
In equation (1.5) 'p' is the total number of colors in the color map. The default 



value for color maps used is p = 256, and the default color is Gray Scale. The user 
of the software can choose from a list of color maps, or can create his own color 
map. 

1.2.2. Unique Findings from the Experiments. As the software can evaluate almost 



any mathematical expressions to generate the fractals of type equation (1.3), we 
noticed certain interesting findings. 

The most interesting phenomenon we noticed is the existence of 'embedded' 
fractals while using Transcendental expressions, that is fractals, which are cre- 
ated by expressions with transcendental functions in it. The first observation was 
l tan(z) 2 + c' fractal ( figure [8~| on page 



inside it (figure 18 on page 19). 



12 



)contains L z 2 + c' or Mandelbrot fractals 



There are Mandelbrot fractals inside l (tan(z)) +c' fracta 



20 



on page 



20 



19 



on page 



19 



. Figure 

depicts the embedded 



illustrates a (cotan(z)) + c fractal, and figure 
Mandelbrot fractals inside it. 

The second observation was, that the Mandelbrot set is capable of 'Shape Pre- 
serving Transformations', like rotation, scaling, and translation, just by multiply- 
ing the z 2 term by different (complex) factors. We have established a 'generic' way 
to transform polynomial fractals without a z 1 term. 



1.3. Generalization of the Transformation Function. Equation (1.1) can be 
further generalized into:- 



F : z <- z n + c; n G N 

where N is the set of Integers. 

Further generalization can be done by:- 

F : z <- z n + c : fi6K 



1.6) 



;i-t) 



where R is the set of real numbers. The resulting sets are known as "Multibrots" 
family. In particular, [3] and jl] consider these maps in details. 
Figure [T] on page [4] shows a "Multibrot" fractal with n = 4.0. 
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Figure 1. Multibrot, with n 



2. Invariance : Shape Preserving Transforms on Fractal Shapes 

We know that the only shape preserving transforms are:- 

(1) Translation 

(2) Rotation 

(3) Scaling 



In this section we review the integer multibrot map (1.6), under the influence 
of shape preserving transformations. Our observation is that Standard shape pre- 
serving transformations works on the function family of 'Multibrots ' with positive 



'n' (equation (1.6) , n > 0). We note that the family with n G Q does not preserve 



shape under rotation, simply because the transform cease to become a function :- 

z <- + c =>• z <- h \fz°- + c 

The expression b y/z^ is not a function, and therefore no 1-1 relationship between 
the old point, and mapped point can be established. 

2.1. Shape Preserving Transformation Theorems. We show that the point 
which is now current transformed u z y , has the same trajectory of the earlier i z\ 
So, if V was included in the set, so would be with modification in fc's value, 
and, if V was not included in the set, so would not be a z\ That would establish 
the transformation properties. 

Theorem 2.1. A variable transformation of the form z : z — a translates 
a Multibrot Map. 

Translation can be achieved by :- 

T t (z) : z <- (z - a) n + c ; a{x,y) G C ; n > (2.1) 

The result is the set image shifted co-ordinates to a = (x,y). 

One point to be noted here is that as the radii of convergence is k, if the point a 
resides further than k, then we need to increase the radii appropriately, for seeing 
the whole set's image. 
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Proof of the Translation Theorem. To prove this we note that all points which were 



mapped as z according to equation (1.6), are now mapped as z + a in the new 



mapping with equation (2.1 ). That would mean that if z 6 M for equation (1.6) 



then z + a G M' for equation (2.1 ), and vice versa. We take a point *c = c + a and 
hence t z = c + a. 



Clearly at 1st level the original transformation can be written as:- 



Z\ 



c n + c 



(2.2) 



Compare with 1st level transformation using equation (|2.1|):- 

t. 



*z - a)" + 'c 



('c-af + 'c 



c + a - a) +c + a = c+c + a 



where t z and *c signifies the transformed z, c. 
So we have established that:- 

'^i = + a 

Now, let it be true for some I > 1. Then we have:- 



(2.3) 



z t = zi + a 



We have:- 
and 



Zl+l = z? + c 



(2.4) 



Zl+l 

Replacing the value of l z\ we have:- 



zi 



r + 'c 



Vi = ( 



Zi + a - a) n +* c 



[z? + c) + a 



hence:- 



zi+i = z i+ i + a 



(2.5) 



Using equation (2.3) as basis and then having equation (2.5), implies that using 



principle of mathematical induction, we can say all points which were mapped as 



z according to equation (1.6), are now mapped as z + a in the new mapping with 
equation (2.1). That establishes the translational property. □ 



Theorem 2.2. Transformation z : e %e z rotates a Multibrot Map. 

A "clockwise" rotation of angle p = Qjin — 1) can be achieved by:- 

T r (z) : z <- e w z n + c (2.6) 

Proof of the Rotation Theorem. To prove that we note again:- zq = c and t zo = 

ce \ n -i' _ 

We have:- 



c n e { 



?i W — -W— nib 



) + ce {n -^ = e { ^\c n + c) 



Hence, using equation ( |2.2 ) 



2i 



3*- n-l ' 



2l 



(2.7) 
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Assume now that for some I > 1 



Zi = e^^'zi 



Clearly we would have:- 



"zi+i =e lu ("z l )+' "c=e 



/ — iO \ f —iO \ / — i& \ 

^ n-1 ' -|- C6 n_1 = p\n—i'i 



Then, using equation (2.4) we get:- 

Zl+l 



( — ) 



Zl+l 



(2.8) 



Using equation (2.7) as basis and having equation (2.8) implies that, using 



principle of mathematical induction, we can say all points which were mapped as 

. . / — iB \ 

z according to equation (1.6), are now mapped as e {n - 1 'z in the new mapping with 



equation (2.6). That establishes the rotational property. 



□ 

Note that for "Multibrots" with n < this formula holds true. As the term 
p < OVn < , this rotation would be in "anti-clock-wise" direction. For n > 
, this rotation is in "clockwise" direction. Also note that, for "Multibrots" with 
fractional n > 0, rotation leads to altogether different type of images. They are 
not invariant under rotation. 

Theorem 2.3. Transformation z : az scales a Multibrot Map 

A Scaling of factor a can be achieved by :- 

T s (z) : z^az n + c; a = a { ^ } ; a G R + ; n > (2.9) 
The resulting figure is a scaled version of the original Map, by a factor of a. 

Proof of the Scaling Theorem. To show this, we start with zq = c and f zo = ca^ 1 ^^ 
as initial points. 
Clearly then:- 

* Zl = a (V) +' c = a(ca ( ^ y ) n + ca { ^ = a ( ^\c n + c) 



hence, using equation (2.2) 



Z\ = a Kl ~ n 'zi 



(2.10) 



Now, assume it is true for some I > 1 



t (^-) 
zi = a Kl -"'zi 



Clearly we can say:- 



Vi = a(V) +* c = a(a ( T^ } ^) n + ca^ = a { ^\z? + c) 
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Then, using equation (2.4) we get:- 



Zl+l 



(—) 



Zl+1 



(2.11) 



Using equation (2.10) as basis and having equation (2.11) implies that, using 



principle of mathematical induction, we can say all points which were mapped as z 
according to equation (1.6), are now mapped as a^^^^z in the new mapping with 



equation (2.9). That establishes the scaling property. 



□ 



Note that, "Multibrots"with fractional n > are invariant under scaling. 



2.2. Generic Shape Preserving Transformation. The most generic shape 
preserving transformation would be, then 

T(z) : z <- a(z-b) n + c ; a = \a\e id , b(x, y) E C (2.12) 

The transformation of equation ( |2.12[ ) translates the origin by b = (x, y) amount, 
then rotates the map by 0/{n — 1), and then scales down the map by lal*- 1 ^. 

2.2.1. Shape Preserving Transformations For Polynomial Fractals. Polynomials 
are important due to Taylors Series expansion available for any infinitely differen- 
tiable function. The generic formula for a polynomial transform is:- 



T(z) 



cijzi + c 

3=0 



P(z) + c 



(2.13) 



We note that while equation (2.1) can be immediately applied to a polynomial, 



due to variations in power transforms like rotation: equation (2.6) and scaling: 



equation (2.9) can not be used on equation (2.13). Hence, only with the translation 



part, we can redefine equation (2.12) for polynomials as 



T(z) : z <- P(z - b) + c ; b(x, y) E C 



(2.14) 



2.2.2. Rotation and Scaling of Truncated Polynomial Fractals. We can not consider 



polynomial with a linear z 1 term because then the results from rotation (2.6) and 
scaling (2.9) becomes undefined with n = 1. 



Definition 2.1. Left Truncated Polynomial. 

A Left Truncated polynomial is defined as:- 



P(z,l) = J2^z k ; l>2 



(2.15) 



k=l 



We define a truncated polynomial is a polynomial with at least no linear term, 
(and possibly constant term) to be defined formally in a later section, with equation 
( |2.15[ ). We can drop the requirement of not having a constant term, because in 
that case z° , and the equations (2.6) and (2.9) remains defined. 
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We have found out from equations (2.6) and (2.9) that they affect different 
powers of 'z' differently. But assume, we still want to find a way to rotate, and 
scale a polynomial fractal image appropriately. For notational convenience, we 



treat P(z) as a vector, with components Xj 



3 ■ 



3" ' 



P = [a jZ j ] 

We note that, to generate an uniform scaling, <7j ^ Oj. 

s = N 

We note that, to generate an uniform rotation, pi ^ pj. 



(2.16) 



r = [ Pj ] (2.17) 
It is easy to see that the for the 'effect' of them all should be equal to some 



predefined value V 
rewriting, we get:- 



u 



°3 



u 



1-3 



In the same way, to get a fixed angle of rotation '#' :- 

Pj = e^ 9 

Hence, a properly scaled and rotated original Polynomial fractal P+c transform, 
with translation, would have components:- 



"*p = [ T x 3 ] = [« l -V 



(2.18) 

Using equation (2.18) we can modify and shift any 'non- infinite' (j ^ oo ) 
truncated (j ^ 1 ) polynomial fractal, such a way that the original shape is 
preserved. So, the set image remains 'invariant' under equation ( 2.18[ ). 

Note that equation (2.18) can handle n < 0, but not the scaling part, as 'fc' 
needed to be adjusted. 



2.2.3. Demonstration of Invariance. We take the equation of : 

T(z) 



z 2 + z 3 + c 



(2.19) 



The image of equation (2.19) can be seen in figure |9jon page 16 The image of 
this can be rotated clockwise ti/2 radian, by multiplying the \2 by i and xs by 
e tn . Hence, 



l T(z) 



%\z 



r 



[z-lf + c 



(2.20) 
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Equation (2.20) is a clockwise n/2 rotated, and right shifted by 1, 



equation (2.19), which is clearly established by the figure 12 on page 17 



image of the 
To further 



establish the result, we present the second equation:- 



T(z) 



e^Hz 



lf + i{z- if + c 



Equation (2.21) is a clockwise tc/A rotated, and right shifted by 1, 



(2.21) 
image of the 



equation (2.19), which is clearly established by the figure 13 on page 17 



2.3. Self Similarity. With the previous set of transforms in mind, we focus on 
what is observed self similarity in the Multibrot fractals. 

The tool we seek is the magnify a small section of the fractal shape, and finding 
how it should look. We note that this has to deal with the transform scaling ( 
theorem 2.3 ). But that won't be all. 



The point with respect to which we need to magnify or zoom first needs to be 
made the origin of the co-ordinate, and then magnified , or scaled. That would 
make the transform looks like translation followed by scaling. But that is not so. 

We established that the direction of the co-ordinate system matters for fractals, 
as in Multibrot with fractional powers does not show remain same under an induced 
co-ordinate rotation. Therefore, the transform should also take care of aligning 
the co ordinates properly, via rotation. 

Definition 2.2. The Zoom Transform. 

Zoom of a fractal image can be done using :- 

z <- se^f(z -a) + c (2.22) 

where s is a positive constant factor, a = (x n ,y n ) a complex number to have the 
co-ordinate origin shifted to, and <fi denotes the rotation angle. 

Conjecture 2.1. Existence of Self Similar Embedded Fractals. 

A fractal of the type :- 

z <- f(z) + c 

has quasi self similar embedded fractals under magnification iff it is invariant under 



the equation (2.22) or the zoom transform. 



3. Generic Analytic Family and Embedded Fractals 

It is to be noted here that all geometric transcendental functions can be ex- 
panded in Taylor series, and hence a transform of the form: 

T : z <- G(z) + c 

where G(z) is a geometric transcendental function, qualifies as an infinite poly- 
nomial map. 



10 



NABARUN MONDAL AND PARTHA P. GHOSH 



Figure 2. Fractal of cos(z) — 1 + c , k is about 8 

i 

Figure 3. Fractal of cos(z) — 1 + c showing Inner Mandelbrot, k 
is about 8 



Figure 4. Fractal of sin(z 2 ) + c , k is about 8 



We have noticed an important property of such maps, that they have embedded 
Multibrot fractals inside small regions within them. Some of them we have already 
shared in the introduction section. We furnish some more examples here:- 

The fractal of l cos(z) — 1 + c' (figure [2] at page 10 ) contains Mandelbrot 



fractals inside (figure [3] at page 10 



(2) The fractal of 'sm(z^+ c' (figure |4j at page 10) has Mandelbrots inside 
(figure [5] at page [IT] ) , and looks like a deformed eaten away Mandelbrot 
image overall. 



PROPERTIES OF POLYNOMIAL AND TRANSCENDENTAL FRACTALS 



11 




Figure 5. 
about 8 



Fractal of sin(z 2 ) + c showing Inner Mandelbrot, k is 



Figure 6. Fractal of sin(z 4 ) + c , k is about 8 



J* 



I 



I 



Figure 7. Fractal of sin(z 2 )+c showing Deformed Inner Multibrot 
with n=4 , k is about 8 



4 ) + c' is (figure |6] at page 11) virtually same 
4 (figure [l] at page [4]) , and has Multibrots n = 4 



(3) Increasing the power, \s///(.; 
as a Multibrot with n 
inside it (figure [7] at page [TT 

(4) What more, the fractal of '6(z — sin(z)) + c' shows up just like a Multibrot 
of n 



3 (figure [8] at page 12). 

These observations lead us to think that for those maps, the original transform 
function can probably be approximated by the Multibrot type function. This is 
the topic of discussion in this section. 
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* 



Figure 8. Fractal of 6(sin(z) — z) + c Resembles a Multibrot with 
n=3 , k is about 8 



3.1. Asymptotic Bounds : Bachmann Landau notation. What we infor- 
mally seek is that in the neighborhood of a point, can we replace a function by 
another much simpler function, because, they behave approximately the same way. 
This concept is already available as asymptotic bounds proposed by Bachmann [6] , 
and reintroduced by Landau [7j, for real functions, with a parameter x. 

Definition 3.1. Big Theta : 6 . 

Let there exist functions f(x),g(x) and constants k\ > 0, k 2 > such that :- 

hf{x) < g(x) < k 2 f(x) ; x x crit 

then, 

g{x) G e(f(x)) (3.1) 



The equation (3.1) is the Bachmann-Landau '0' notation, mostly used in com- 
puter science (Knuth [8]) for getting approximate bounds for functions (runtime 
analysis of algorithms). With simple modification, complex analytic functions can 
be included in 6 notation. 

Definition 3.2. Complex Generalization of Big Theta. 

Iff for F(z), G(z) and constants k\,k 2 G (0, k) at neighborhood of a point z c G C:- 

h\F(z)\<\G(z)\<k 2 \F(z)\; z -> z c (3.2) 

then, we define:- 

G(z c ) G 0(F(O) (3.3) 
As both F(z), G(z) are analytic, then, there are infinite points in the neighbor- 



hood of z c where equation (3.2) will be true. On that neighborhood, 1-^(^)1 will 



approximate |G(^)| with arbitrary accuracy. 

Theorem 3.1. There is a relationship for Complex Polynomials. 

Let 

U 

f(z) = ^a k z k ■ a k eR 

k=l 

Then, if this holds:- 

lim ]— — = 

z^tz c \z\ m 
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then 

lim f(z) G 0{z m ) 

z—>z c 

In particular v- 

y\ z \ -> oo =^ /(z) G 0(,2 U ) 

and, 

V|z| -> /(z) G e(z l ) 



Proof of the theorem \3.1\ . We rewrite the polynomial as:- 

v fc 



f(z) = z m J2^ 



\z 



k 



k=l 

Taking mods we get:- 

l/(*)l = N m lE fl *tiWi 

which immediately gives the result. □ 

Theorem 3.2. The relation is true over iterations over polynomials. 

Let f(z) be a polynomial function with 

lim f(z) G e{g(z))i 

where g(z) = z m Then, 

lim p{z) G ®(g n (z)) 



Proof of the theorem 3J2 . We rewrite the polynomial as:- 

u u 

f(z) = z m J2 a * Z 

and then the f n (z) becomes:- 

P(z)=z mn J2 b *- 



z u 

k=l 



un u 
Z K 



Z' 
k=ln 



Taking mod, the result follows immediately. 



□ 
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3.2. The Principle of Dominance. What we informally seek, is when we "zoom 
in" to a polynomial fractal, what "inner" fractal type we may find. If we find a 
fractal of type f c {z) = f(z) + c at a point l Zd inside a transformation function 
T c (z) = T(z) + c, then we say that the function f c (z) must have "dominated" the 
function T c (z) at point z&. 

Theorem 3.3. Relationship © can be used to find escape time fractals. 

Let the transform function z <— T(z) be:- 

T = f(z) + c 

f(z) being a polynomial with lim^^ f(z) G 0(g) . Let then:- 



and 

so that:- 



S={z: lim \f n (z)\ <oo} 

n^ooz^z^ 



S e = {z : lim \g n {z)\ < oo} 



S — *S*( 



e 



Proof of the theorem 3.3 We note that from theorem (3.2) :- 



This would imply that:- 
which would imply that:- 
also 

which would imply that:- 
combining we get Sq = S. 



lim f n E Q(g r - 



in < km 

SCS e 
\r\>Ki\g n 
Se^S 



□ 



Theorem 3.4. Dominance : Relation f G 0(g) guarantees a 'sub- fractal'. 

Lets lim 2 ^ d / G 0(g). Then, there for the escape time fractal generated from 
the transform function:- 

z <- f(z) + c 
there will be a embedded sub fractal of the form:- 



In general for f \z) = P(z,l) , a truncated polynomial (definition 2.1), there would 
be sub fractals of type z l + c. 
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Theorem |3.3| makes proof of this theorem a triviality. However, another proof 
can be furnished from the perspective of scaling transform. We demonstrate the 
proof here. 



Proof of the theorem \3Jh ■ We assume that:- 

l\g(z)\ < 1/0)1 < u\g(z)\ with u, I > 
If that is the case, then we know that the fractal:- 

z lg(z) + c 
would be bigger in size than the fractal:- 

z ug(z) + c 

if the fractal z g(z) + c exists due to the scaling theorem (2.9). 

That would mean that the fractal z <— f(z) + c is to be bound between the 
two fractals of the same shape, which is only possible, iff the shape of the fractal 
z f(z) + c is similar to that of z g(z) + c, which proves the theorem. 

□ 



3.3. Implications of Dominance. The theorem (3.4) suggests, that there will 
be Multibrots found in the vicinity of truncated polynomials. In this section we 
elaborate what observations we have made, and how those observations fit into the 
property of dominance. We demonstrate simple truncated polynomial fractals. 

3.3.1. Demonstration of Dominance. To illustrate the result, we exhibit two trans- 
formation functions. The first transformation function:- 

T 2>3 (z) :z 2 + z 3 + c 



has Mandelbrot [z + c) fractals inside it. Figure |9j on page 16 shows the T 2 ^(z) 
and inner (z 2 + c) in the next image. 
The seccond transformation function:- 



24,7,-io(-z) : z 4 + z 7 -z 10 + c 



has Multibrot (z A + c) fractals inside it. Figure 14 on page 17 shows the T 4)7 _ w (z), 
and inner (z 4 + c) in the next image. 

These two observations shows validity of the dominance property. 

3.3.2. Implication for Functions Expandable with Taylor Series. The most impor- 
tant result of Dominance lies with the Taylor Series expansion. For a function 
f(z) which is capable of expansion by Taylor Series, "multibrot" type fractals can 
be generated by removing lower powers, or multiplying or dividing the function by 



z k . 



We note down the expansion of sin(z) 



z 3 z 5 



sin(z) = z H — - — ... (3.4) 

v ; 3! 5! K J 
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Figure 9. The Polynomial Fractal of z 2 + z 3 + c 



4 



Figure 10. z 2 + z 3 + c zoomed, showing existence of z 2 + c fractal inside 



Figure 11. z 2 + z 3 + c, zoomed, showing clearly + c fractals 
and the expansion of cos(z) :- 



z 2 z A 



cos(z) = l- - + --... (3.5) 
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Figure 12. Image of (z — l) 2 — (z — l) 3 +c, showing the invariance 
under translation and rotation. Cross reference with Figure [9] on 
page [16 



# 



Figure 13. Image of e^\z — l) 2 + i(z — l) 3 + c, Cross reference 



with Figure 12 on page 17 




Figure 14. The Polynomial Fractal of z 4 + z 7 — z 10 + c, notice the 
similarity with figure : [T] on page [4] 



In this subsection for brevity we have used the informal G ~ F to denote 

g g e,(F). 

Now clearly then due to dominance effect:- 
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Figure 15. The Polynomial Fractal of z 4 + z 7 — z 10 + c, zoomed, 
showing how the borders look like a Multibrot with n = 4 

Figure 16. ,2 4 + z 7 — z 10 + c zoomed, showing z 4 + c fractals inside, 
compare with figure : [T] on page [|] 

(1) 

cos(2i) — 1 + c ~ — + c— >■ z 2 + c 

(2) 

6(sm(z) — z + c) ~ z 3 + c 

(3) 

sin(z n ) + c~ z n + c 

(4) 

z^mtz") + c ~ z n+k + c ; fceK 

(5) And finally we explain, the initial tan(z) fractals, which served as the 
starting point of this paper, in the next subsection. 
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Figure 17. (tan(z)) 2 + c Fractal 




Figure 18. (tan(z)) 2 + c zoomed, showing Mandelbrots inside 




Figure 19. (cotan(z)) 2 + c Fractal 



3.4. The tan(z) Fractals. Figure 17 on page 

111 

shows a (tan(z)) + c fractal, 
and figure fLS] on page [19] shows embedded Mandelbrots inside it. 



hence, 



z 3 2z 5 

tan(z) = z H 1 h ... 

V ; 3 15 



(tan(z)f + c ~ z n + c 



(3.6) 



Now we proceed to explain the observation of '(tan(z)) 2 + c\ Figure 



19 



shows a (cotan(z)) + c fractal, and figure 



page 

Mandelbrots inside it. 



20 



on page 



20 



on 



shows embedded 
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Figure 20. (cotan(z)) 2 + c zoomed, showing Mandelbrots inside 




Figure 21. T g (l, 3, 2) Fractal, note the similarity with figure 19 on page 19 



(3.7) 



1 . , 1 z z 3 2z 5 

— = cotan(z) = ... 

tan(z) w z 3 45 945 

We note down that equation (3.7) contains power of z" 1 , and hence, we can not 



directly use the theorem of existence of embedded Multibrots. 
Let us define :- 



g{a,b,z) 



; a, b G C 



and 



(3.8) 



T g (a, b,n) : z <— (g(a, b, z)) n + c 

A fractal of type T g (l, 3, 2) is shown in figure 19 on page 
We state without proof:- 
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(3.9) 



(cotan( y z)) n + c ~ T g (l, 3, n) 



\z\ < 1 



and 



(cotan(z)) n + c ~ z n + c ; |^| > 1 

Theorem 3.5. Existence of Multibrot inside Fractal of T g (a,b,n). 
Fractal of equation family T g (a, b, n) (3.9) would contain Multibrots inside. 
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Figure 22. T g (2, 5, 3) Fractal, 1000 times zoomed, Shows a Multi- 
brot inside with n =3 



Proof of theorem (3.5). We note that when z — > 0, the equation (3.8) gets un- 
bounded. 



Assume k is large enough so that:- 



k > a 



1 



Now, when z — > 1 but z ^ 1 and |z| < 1, we have:- 

z = l-e; e G C ; e^O; e^O 

In that case:- 



lim g(a, b, z) = lim( 

z— >1 ' e— >o 1 — e 



1 



-) 



That approximates into :- 

lim(a(l + e + e 2 + e 3 ...) - 
Ignoring the higher powers of e we get:- 



a(l + e) 



1 - e 



(3.10) 



g(a, b, z) = 2a — z(a + -) ; z — > 1 ; 1 



And hence equation ( 3.9[ ) can be written as:- 



T g (a, b,n) : z (^(a, 6, z)) n + c ~ (2a — z(a + -)) + c 



(3.11) 



(3.12) 



But, equation (3.12) is a 'shape preserving transform' of equation (2.12). Hence, 
the fractal of equation (3.12) is expected to have Multibrot fractals inside. □ 



A zoom of fractal of type T g (2, 5, 3) is shown in figure 22 on page 21, revealing 
the Multibrot with n = 3 inside it. 

Hence we complete the explanations of the observations what we have stated in 
the beginning of the current section. 
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4. Summary and Future Work and Thanks 

This paper started as a celebration of what Mandelbrot's work has been. We 
started working on creating a software tool to produce 'only' Mandelbrot set in 
color, and ended up creating a software tool capable of evaluating any Mathemati- 
cal expression, to generate escape time fractals, in any color. While experimenting 
with this tool, we came up with heather to unknown, unexplained phenomenon. 
The phenomenon of shape preserving transform, and the phenomenon of embed- 
ded inner fractal images, were so fascinating that we tried to explain them, so as 
to generate a preliminary theoretical basis. The usage of transformation vector to 
preserve shape and transform any finite truncated polynomial fractal is proven. 
The dominance phenomenon we borrowed, and reused from Computer Science, 
and shown that it is a very powerful mechanism to predict existence of embedded 
inner fractals. While dominance principle shows promise, the principle is fairly 
new at least in the fractal geometry domain. Further theoretical research work is 
needed to take the preliminary concepts discussed in this paper, into mainstream 
fractal geometry. 

We would specially like to thank Carl Johansen, whose open sourced Mandelbrot 
viewer we borrowed, without which, not a single line of this paper could have been 
written. 
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